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In theories with flat directions containing vortices, such as supersymmetric QED, there is a vacuum 
selection effect in the allowed asymptotic configurations. We explain the role played by gauge fields 
in this effect and give a simple criterion for determining what vacua will be chosen, namely those 
that minimise the vector mass. We then consider the effect of vacuum selection on stable (BPS) 
non-topological vortices in a simple Abelian model with N = 2 supersymmetry which occurs as a low 
energy limit of Calabi-Yau compactifications of type II superstrings. In this case the magnetic flux 
spreads over an arbitrarily large area. We discuss the implications for cosmology and for superstring 
. inspired magnetic confinement scenarios. 
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I. INTRODUCTION 



Flat directions of the scalar potential are a feature of a large class of supersymmetric (SUSY) theories, for instance 
in Abelian theories where the gauge symmetry is broken with a Fayet-Iliopoulos (FI) D-term. In the last few years 
these have undergone intense study, in particular in the context of effective actions for SUSY non- Abelian gauge 
| theories and confinement (see e.g. [|, §))■ 

. A generic consequence of SUSY theories with a D-term is the formation of cosmic strings |j, Q]. However, the string 
' solution is richer than the usual Nielsen-Olesen (NO) vortex |s| ^| for two distinct reasons. Firstly, the existence of 
ON , flat directions means that there is a moduli of degenerate vacua and, at first sight, it is unclear which of the vacuum 
■ states give rise to stable cosmic string solutions. Secondly, the presence of supersymmetry gives rise to fermion zero 
modes in the string core [Q, rendering the string current-carrying. This typically leads to the formation of vortons, 
long-lived remnants whose observational signatures can severely constrain, and in some cases rule out, such models 0. 

Topological strings arising from the bosonic sector of an N — 1 SUSY model have been analysed in ]3| , showing that 
even though there is a one-parameter choice of vacuum for the moduli field, only one choice permits cosmic string 
solutions which could be stable. Their conclusion was based on the behaviour of the modulus field well outside the 
core of the vortex. This vacuum selection effect is generic in Abelian theories with flat directions in that, as we will 
show, it is driven by the tendency to minimise the vector mass at the core of the vortex. Thus, one expects it to affect 
(— ( non-topological as well as topological strings, and it could conceivably improve their stability 

In this paper we analyse this vacuum selection effect in the context of (BPS) non-topological strings. In particular, 
we investigate the nature of BPS vortices in N = 2 supersymmetric QED with two hypermultiplets of opposite charge. 
This model has been analysed previously in connection with Calabi-Yau compactifications of type II superstrings |§ , 
where the existence of vortices is linked to the confinement of magnetic charges. We will find that, due to the vacuum 
selection effect, the structure of the BPS vortices is identical to that of so-called semilocal strings^). 

Apart from its physical significance, this model has the advantage that all calculations can be performed explicitly, 
which is important because the stability of semilocal strings is somewhat counter-intuitive. We should stress that 
although the vacuum manifold is simply connected, the gauge boson is massive and magnetic flux is topologically 
conserved and quantized; the vortices analysed here are bona-fide BPS states carrying one unit of magnetic flux. 

And yet we will find that, in practice, there are no stable vortices in this model. To be more precise, the vortices 
are only neutrally stable (which is consistent with the BPS condition) and degenerate with a whole family of other 
BPS "magnetic flux tubes" of arbitrarily large width. The slightest perturbation excites the zero mode and makes the 
vortex expand. The vacuum selection effect does not single out the narrow Nielsen-Olesen vortex over the other, more 
extended, BPS flux tubes. In a cosmological context, this suppresses vorton production. In the context of superstring 
compactifications, it casts serious doubts on the mechanism of magnetic charge confinement proposed in 

The paper is organized as follows. We first review the vacuum selection effect in the model of ||. We explicitly 
show that the resulting cosmic string saturates the Bogomol'nyi bound, so is strictly stable, and that it is no accident 
that the selected vacuum is the one minimising the vector boson mass. We argue that this gives a generic criterion 
to determine which vacuum or vacua are selected in a given theory, provided all charges are equal in absolute value. 

We then extend this idea to N = 2 QED with two hypermultiplets of opposite charge, using a Bogomol'nyi argument 
on the corresponding solution pToj ] . We again show that it saturates a Bogomol'nyi bound, so is stable, but the presence 
of a cylindrically symmetric zero mode changes the dynamics completely and in particular the string expands to an 
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arbitrarily large area in finite time [O, [12| |l3j . This occurs even though the vector mass is non-zero and the magnetic 
flux at infinity remains quantised, setting a constraint on viable models of confinement of magnetic charge in these 
N=2 Abelian theories. Finally we discuss the effect of SUSY breaking on the solutions and speculate on the resulting 
cosmology. 



II. VACUUM SELECTION EFFECT AND MINIMUM GAUGE MASS 

We begin by reviewing the vacuum selection effect described in ||. Penin, Rubakov, Tinyakov and Troitsky 
considered a four-dimensional model with two N — 1 chiral superfields of opposite charges coupled to a U(l) vector 
multiplet, and a FI D-term which gives the symmetry breaking potential. The action of the bosonic sector is 

S =fd 4 x [\D^ + \ 2 + \D^\ 2 - - V{4>+A-j\ 

where is a U(l) gauge field and <f>± are complex scalars of opposite charges (q± = ±q), 

Dfj,(f>± = (d f _ l ±iqA fl )(j)±, F M „ = d^A v - d^A^. 

The solutions for the vacuum manifold (V = 0) are 

|0+| = rj cosh u = v + , \4>—\ = r)smhu = V— , (2) 

where u parametrises the moduli space. After symmetry breaking the physical spectrum consists of two massless 
scalars, the Goldstone boson and the modulus field; a massive scalar particle m 2 = 2A?7 2 cosh2w; and a massive 
vector particle m 2 — 2q 2 r/ 2 cosh2w. Their masses depend on the choice of vacuum. The minimum values fh s and fh v 
correspond to u = (<f>- — 0), where we recover the Abelian Higgs model. This has static vortex solutions, the NO 
strings §. 

We want straight, static vortices along the z-direction so we drop the t- and z-dependence and set A t — A z = 0, 
defining B = d\A^ — 82 A\. We can set r\ = q = 1 by rescaling, and we consider the Bogornol'nyi limit A = 1 this 
effectively makes the model N=2 supersymmetric with one hypermultiplet, see also Jl^, |l4|, Kjfo. The conditions for 
finite energy include D^<f> + — > 0, D^cj)- — > 0, B ~ > faster than 1/r and this introduces a correlation in the phases of 
the fields at infinity: fields with charge q a wind as e mqa9 and the gauge field tends to a constant, Ae — ► — n, leading 
to the quantisation of magnetic flux, J d 2 B = —2wn. 

In principle, one can try to construct vortices tending to any of the vacua (|J) as r — > 00. However, in j^] it was 
shown that the only static solutions are those with u = 0. Any other choice of boundary conditions leads to what is 
effectively an unstable vortex that tends to this one. The instability is very mild, and the surviving string is a NO 
string. We will now show that this string is stable. 

In what follows we consider cylindrically symmetric configurations with n = 1 : 

0+ = f + (r)e w 0_ = ./L(r)e- l V A Agd0 = a{r)d6 

f±(r) are real functions, with /+(0) = /-(0) = a(0) = and, /±(oo) = v±, a(oo) = — 1. 7--dependent phases e 1 ^^^ 
for <f>± are possible but minimum energy requires d r ^p± = 0, so we will ignore them. A is a real constant. 
To prove stability we obtain the Bogomol'nyi equations 

/iT— / ± =0, ^- { fl-f_-l) = Q. (3) 
r r 

Note that, if /_(r) = 0, eqns. @ are the standard Bogomol'nyi equations for the Abelian Higgs model. 

It can be seen that the only solution for /_ satisfying the boundary conditions at the core is /_(r) = 0; thus 
the NO solution saturates the Bogomol'nyi bound automatically. This shows that the string is stable, since it is a 
global minimum of the energy, and excludes the possibility of cylindrically symmetric zero modes. In fact it turns out 
that the restriction to cylindrical symmetry is unnecessary: always vanishes on the solutions to the Bogomol'nyi 
equations (see, e.g. JpJ). The only remaining solutions are the standard NO vortices. 

As explained in pL the key point in understanding the vacuum selection effect is the effective separation between 
the dynamics of the magnetic core and of the modulus field outside the core (this separation is only true for vortices, 
as it is a consequence of the properties of massless fields in two dimensions, the transverse dimensions of the vortex). 
Far from the core the magnetic field vanishes, there is no potential energy, and the modulus field (properly normalised) 
becomes massless, with solutions logarithmic in r - unless u — 0. Thus u = is selected. 
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We can see this in another way: well outside the core the scalar fields can move along the moduli space from u = 
to their asymptotic values, with no appreciable cost in energy. Indeed, outside the core, 

E~ J d 2 x{D^(j) + ) 2 + (D tl (j)-) 2 ~ J d 2 x{d r u) 2 cosh2u 

and the minimum energy configuration u(r) that interpolates between any two distinct vacua, u(R\) — u\ and 
u{R*t) = U2, if i?2 >> r core , is of the form 

rU2 "1 2 

I duvcosh2u 

mi 

which can be made arbitrarily small as R 2 — > 00. Thus, the dynamics at the core are decoupled from the behaviour 
of the moduli fields; the core of the vortex is effectively free to choose its boundary conditions in order to minimise 
its own energy. At the core (where the magnetic field is non-zero) the field 0_ is suppressed and the core becomes 
identical to that of a NO string. But the previous argument also shows that the minimum energy is unattainable 
unless I — (u\ — 1*2), and this is a condition for static solutions. Thus, the vacuum selection actually occurs at the 
core. 

The authors of || offered no insight into why a particular vacuum is selected, but it is clear that 0_ = is chosen 
because it minimises the mass of the gauge field, m v . This is an important consideration at the core of the vortex, 
where the magnetic field is concentrated in a region of order m~ l . Magnetic field lines repel, so a lower m v means a 
larger magnetic core. Since the total magnetic flux is quantised, this lowers the energy at the core and, therefore, the 
total energy. 
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III. BPS VORTICES IN N — 2 SQED AND VACUUM SELECTION 

We now consider two N=2 hypermultiplets h a , a = 1, 2 with charges q a = +q, —q coupled to an N=2 Abelian vector 
multiplet [[l6| . This model has been analysed before as a toy model for the low energy effective action of Calabi-Yau 
compactifications of type II superstrings by Greene, Morrison and Vafa[||, who identified magnetic configurations due 
to the wrapping of D-branes on cycles. In the low energy theory these would appear as monopole-antimonopole pairs 
joined by non-topological magnetic vortices providing magnetic confinement. 

In [h6 it was shown that, in the absence of FI terms, such vortices were unstable. In an attempt to stabilise the 
vortices, here we consider the addition of FI D-terms, the only modification still compatible with N=2 supersymmetry 
of the effective action. We will be partially successful in that the instability turns into a zero mode, but as we will 
show this is still not enough to provide confinement. 

The model has a global SU(2) symmetry that rotates the two scalar fields in each multiplet, and we introduce 
a triplet of FI D-tcrms k ■ D. In this case the answer can be more involved: since the vacuum manifold is simply 
connected there is no unique vacuum that minimises the vector mass and it is not clear that the selected vacua can 
lead to vortices. We take k = (0, 0, i |g|u; 2 ) , without loss of generality; note that this breaks the SU(2) symmetry. 

Eliminating the auxiliary fields and setting to zero those fields that have zero vacuum expectation values, simplifies 
the model considerably. With a suitable rescaling hij — > u>hij, x^ — > x^/tuq, — > tuA^, the energy of static, straight 
vortex configurations can be written as 



2E 



l^iil 2 + \D^h 12 \ 2 + \D^h 21 \ 2 + \D^h 22 \ 2 + X -B 2 



d 2 x 

{H x l - H 2 + 1/2) 2 + (H, 1 + H 2 ) 2 + (i Hi - i H 2 ) 2 } (4) 



where 0^ = 0^ + i(q a /q)A^ and Hj 1 = -(q a /2q) h* al h aj . 

There is another SU(2) symmetry between (hn and h 22 ) and (h\ 2 and h 21 ) which is preserved after adding the 
D-term. For finite energy we require D^h ai = as r — > 00, which correlates the phases of the multiplets at spatial 
infinity, giving h a ~ e m i^d an( j ^he quantisation of magnetic flux from Agd0 ~ —nd9; the scalars must lie in the 
vacuum manifold (H 2 l = H? = 0, H 2 2 - = \). 

Parametrising the scalars as h ai = r a i e 1 Xai , gives 

Xu ~ X12 = X21 - X22 + 2tott , (5) 
rn 7*12 - r 21 r 22 = 0, (6) 
(rn) 2 + (r 22 ) 2 -(r 12 ) 2 -(r 21 ) 2 = l. (7) 



4 



The vacuum manifold is simply connected, but the set of finite energy configurations has non-contractible loops 
labelled by n, the (quantised) magnetic flux. This situation is familiar in the context of semilocal models (see |U| 
for a recent review), where the gauge field couples to two (or more) scalar fields of equal charges, and there is a 
global SU(2) (or larger) symmetry between the scalars. The stability of semilocal vortices depends on the masses of 
the scalar and the vector particles ]TT| . If m s < m v the strings are stable. In the Bogomol'nyi limit m s — m v , the 
resulting n — 1 "vortices" have a richer structure including a scalar condensate at the core which causes the magnetic 
flux to spread over an arbitrarily large area. The thinnest vortex is a NO string but it is not protected against the zero 
mode that generates the scalar condensate; as a result its energy density becomes arbitrarily close to the vacuum in 
a finite time and the string dissolves |Q. If m s > m v the zero mode becomes an instability and there are no strings. 

Let us consider the vacuum selection effect in model (Q) . The mass of the gauge field is r\ x + r\ 2 + r 21 + r 22 , and 
the minimization of this mass, subject to eqns (^|J7|), predicts h\ 2 = h 2 \ = 0. To prove that this is the case, we rewrite 
the energy as 



2E = J d 2 x[\{D 1 +iD 2 )h ll \ 2 + \{D l -iD 2 )h 12 \ 
+ \(Di + *D 2 )h 21 \ 2 + - iD 2 )h 22 \ 2 
+ [B + {H 1 1 -H 2 2 + l -)] 2 + (H^+H, 2 ) 2 



+ {iH 2 l -iH^) 2 }- j ' d 2 xB , (8) 

and obtain the Bogomol'nyi equations: 

{Dx+iD 2 )hn =0 {D 1 -iD 2 )h 22 =Q (9) 

(Di - iD 2 )h 12 = {Di + iD 2 )h 21 = (10) 

H 1 2 = (11) 

B+ [iV -H 2 2 + l/2] = . (12) 

(Note that setting H x 2 = implies D+R x 2 = which holds if D + hn = D + h* 12 = D + h 2 i = D + h 22 = so the 
choice of ± signs in the gradient terms is not arbitrary). 

The most general configuration for the scalar hypermultiplets compatible with cylindrical symmetry is |l6| ] 

A e = a{r) A r = . (13) 

To analyse the Bogomol'nyi equations we simply note that the pairs of fields (/in, h 2 i) and {h 22 , h* l2 ) behave exactly 
like <j)-) in the N=l model of the previous section: they have charges ±1 and they appear in the square bracket 
of eq. ( |l2| ) with signs +1, -1 respectively. The only way to satisfy ( |l0|) is to have h\ 2 = h 2 \ = 0, which satisfies (|Tl| ) 
automatically. These configurations minimise the vector mass, so we conclude that the "vacuum selection" effect also 
works here (note that by setting the second hypermultiplet to zero the vacuum selection effect gives a topological 
vortex in hn |^5[). 

This leaves eqns. (g), ( |l2|), which are precisely the Bogomol'nyi equations of semilocal strings ^ in (hn,h 22 ). 
These have been studied in ]11|, [l7[ |, so we simply state the results. For n = 1, the non-zero components of the most 
general cylindrically symmetric solution are 

(hn , h* 22 ) = {f{r)e ie , p(r)e lA ) U 

where U is an arbitrary global SU(2) rotation, A is a constant and f(r) and p(r) are real functions with f(0) = 0, 
/(oo) = 1 and p(0) ^ 0, p(oo) = 0. They are related by p(r) = af(r)/r, where a is a real positive parameter which 
labels the solutions. For each a the gauge field is obtained from the condition (£>i + iD 2 )h\\ = 0. 

The lack of winding in p results in the scalar condensate at the core, p(0) ^ 0. Note also that the SU(2) transforma- 
tions between tin and h 22 are not related to the SU(2) transformations coming from the N=2 supersymmetry, since 
they mix elements of the two hypermultiplets. The zero mode is not related to either of these global symmetries, as 
can be seen by the fact that it does not alter the boundaryconditions at infinity. To understand its effect, we consider 



the asymptotic behaviour of the fields far from the core |12j, 18 



f(r) ~ l-i(a 2 /r 2 ) + a 2 (3/8 a 2 -2)/r 4 + ... 

a(r) l + (a 2 /r 2 ) - a 2 (a 2 - 8)/r 4 + . . . (14) 
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There is a one-parameter set of magnetic "vortices" with varying widths. They all saturate the Bogomol'nyi bound, 
so they are all degenerate in energy, but they have different structures: most notably, they have a scalar condensate 
at the core and the magnetic field fall-off is a power of r (in NO vortices the fall-off is exponential). These "vortices" 
can be thought of as hybrids between the NO vortex and the CP 1 lumps, and a can be seen as the width above that 
of a NO vortex: a = corresponds to the NO string modulo SU(2) rotations; but as a — ► oo p(Q) — > 1 and the core 
expands. 

Finally, the zero mode corresponds to a flat direction in the potential for which the vector mass is always at a 
minimum, so the vacuum selection effect does not single out the NO solution over the more extended vortices. 

We conclude that magnetic flux in this model is still not confined to tubes of a definite size. Semilocal strings 
which saturate the Bogomol'nyi bound are at the limiting case between stability (for m s < m v ) and instability (for 
m s > m v ) |jnj [l|. However, |L2| has shown that, in fact, the flux is unconfined and eventually spreads out to 
flux tubes of greater and greater radius. Of crucial importance to the resulting cosmology is the time-scale for this 
relaxation, which is outside the scope of this investigation, but it is clear that vorton production will be suppressed. 

Our conclusions are easily extended to any number of hypermultiplets, as long as their charges are equal in absolute 
value. They also apply in the case of several gauge fields, such as in the low energy limit of type II superstrings 
compactified on Calabi-Yau manifolds, which effectively contains 15 copies of the model analysed here [|J. 

The situation changes when considering hypermultiplets with different charges. In that case vacuum selection may 
be frustrated and the selected vacuum or vacua may not be the ones that minimise the vector mass. The structure of 
the resulting vortices can be totally different from what we described here, with the possibility of binary or multiple 
cores and other exotic effects and it remains a very interesting open question. 



IV. DISCUSSION 



We have investigated cosmic string solutions in SUSY QED with flat directions where the U(l) symmetry is broken 
by a FI D-term. We have argued that, if all matter fields have the same absolute value of the charge, the vacuum 
selection effect acts to minimise the vector boson mass. This includes the original model of j|, the bosonic sector 
of N=l SUSY QED with two chiral superfields of opposite charges, and its simplest extension to N=2 with one 
hypermultiplet, considered in (l5). In both cases the vacuum selection effect leads to topological vortices; we have 
shown that both solutions saturate Bogomol'nyi bounds and are consequently stable. 

We then considered extensions to N=2 SUSY QED with several hypermultiplets of equal \q a \- The case of two 
hypermultiplets of opposite charge was solved explicitly, and there is still a vacuum selection effect, but this case gives 
rise to a semilocal string instead. Again the string solution saturates the Bogomol'nyi bound, but is only critically 
stable as there is a zero mode (which is not an SU(2) rotation). The solutions that saturate the Bogomol'nyi bound 
are parametrised by one or more internal degrees of freedom and are distinct from the string solution, in particular 
their core widths can be arbitrarily large. Note that our N — 2 results generalise to Abelian theories with M 
hypermultiplets provided they have equal |g |; one would still find semilocal strings, but the internal symmetry would 
be SU{M) and not SU{2) Q. Magnetic flux would still be unconfined. 

The situation changes when considering hypermultiplets with different charges. In that case vacuum selection may 
be frustrated and the selected vacuum os vacua may not be the ones that minimise the vector mass. The structure of 
the resulting vortices can be totally different from what we describe here, with the possibility of binary or multiple 
cores and other exotic effects and it remains a very interesting open question. 

As a consequence of supersymmetry, there will be fermion zero modes in the string core for both models. The 
zero modes can be found using the results of Q. For the model of H there will be a single zero mode, left or right 
moving depending on the sign of n, which is a combination of the Higgsino and gaugino. Thus the string will be 
chiral jl9), since it has a chiral current. We can then apply the results of JtJ to this model. If there are no further 
phase transitions to destabilise the current, the model will lead to vorton formation, restricting the scale of symmetry 
breaking |?j . As shown here, this is in contrast with the situation in the N = 2 model, where vorton production is not 
expected to occur. 

However, if we completely break SUSY, for instance via mass terms for the hij fields, we still expect a topological 
string with fermion zero modes in the core. The angular momentum of the resulting current can stabilise loops, giving 
rise to vorton production p0|. The model will then be subject to the constraints of (7). 
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